Anomalous magnetoconductance due to weak localization in 2D systems with 
anisotropic scattering: computer simulation 
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The results of computer simulation of particle motion over the plane with randomly distributed 
scatters are presented. They are used to analyse the influence of scattering anisotropy on the 
negative magnetoresistance due to weak localisation. It is shown that the magnetic field dependence 
of magnetoresistance in this case can be described by the well known expression, obtained in the 
diffusion limit for isotropic scatternig, but with the prefactor less than unity and breaking-phase 
length which differs from the true value. 

PACS numbers: 73.20.Fz, 72.20.Dp, 72.f0.-d 



I. INTRODUCTION 

Quantum correction to the conductivity arises from in- 
terference of electron waves scattered along closed trajec- 
tories in opposite directions. An external magnetic field 
applied perpendicular to the 2D layer destroys the in- 
terference and suppresses the quantum correction. This 
results in anomalous negative magnetoresistance, vifhich 
is experimentally observed in many 2D systems. This 
phenomenon is usually described in the framework of 
quasiclassical approximation which is justified under the 
condition kpl ^ 1, where kp is the Fermi wave vector, 
/ is the mean free path. In this case the conductivity 
correction is expressed through the classical probability 
density W for an electron to retuxn to the area of the 
order Aj;-^ around the start pointEl u 
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This expression allows to calculate the conductivity cor- 
rection at an arbitrary magnetic fieldQ. Analitical expres- 
sions are derived only for some specific cases: (i) in the 
diffusion limit, i.e. at i? ^ Btr and / <C Ijaa^ where l^ is 
the phase-breaking length due to inelastic processes, and 
Btr — 'i'o/27rP,|-|<i>o = irfic/e] (ii) in the high-field Hmit, 
i.e. at B > BtM- 

To our knowledge, this problem has been solved only 
for the case of random distribution of scattering centers 
and isotropic scattering. As a rule, these conditions are 
not fulfilled in real semiconductor structures. First of 
all, in semiconductors the scattering by ionized impu- 
rities dominates at low temperatures. This scattering 
is strongly anisotropic, in particular in heterostructure 



with a remote doping layer. Besides, the impurity distri- 
bution is correlated to some extend due to Coulomb re- 
pulsion of the impurity ions at growing temperatures. In 
the present work the role of scattering anisotropy in the 
quantum correction to the conductivity is investigated 
through a computer simulation. 



II. SIMULATION DETAILS 

To find the values of the mean free path I and probabil- 
ity density W in (|^), we have simulated the motion of a 
particle in the 2D plane with scattering centers in it. The 
plane is a 3000 x 3000 lattice. The scatters are randomly 
distributed in the lattice sites. The whole scatters num- 
ber is 2 X 10''. Every scatter covers seven lattice param- 
eters in diameter. The start point of particle motion is 
chosen in a scatter near the centre of the lattice. We sup- 
pose the particle to move linearly between two sequential 
collisions with scatters (Fig. 1^) . The collisions change the 
motion direction according to given angle dependence of 
scattering probability. A trajectory is considered to be 
closed if after a number of collisions n < Umax it passes 
near the start point at the distance less than d/2. Since 
d has to be small enough, we choose d of about the scat- 
ter diameter. We assume that the particle never returns 
to the area of the start point at n > rimax = 1000 or 
when it escapes the lattice. An estimation shows that 
both assumptions introduce the error in Sa less than one 
percent at 7 = l/l^ > 0.01. 

The simulation has been carried out for two different 
scattering mechanisms: isotropic and anisotropic. For 
the isotropic scattering corresponding to a short-range 
scattering potential, the scattering probability does not 
depend on the scattering angle (the curve i in the inset in 



* email: Alexander.Germanenko@usu.ru 



Fig. p. For the anisotropic scattering we have used the 
angle dependence of scattering probabiUty presented by 
the curve a. It is close to that in heterostructures with a 
doped barrier, where impurities are spaced from the 2D 
gas. The curve a corresponds, for example, to the 2D 
structure with impurity density of about 10^^ cm^^ and 
a spacer of 50 A thick. 

To calculate the conductivity correction in a magnetic 
field, we have modified the expression (||) by including 
the magnetic field in it according to a standart proce- 
dure. The final expression for the conductivity correction 
in the magnetic field can be written as: 
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where summation runs over all closed trajectories among 
a total number of trajectories N (TV = 10^ in our calcu- 
lations), b = B/Btr, Go = e^/27r^?i. Si and k stand for 
the area and length of the i-th trajectory, respectively, 
the exponent accounts for the phase breaking. Note that 
in (2) I denotes the mean free path connected with the 
transport relaxation time. 
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FIG. 1. A fragment of the lattice with scattering centers. 

Lines show closed trajectories. In the inset the angle depen- 
dences of scattering probability for the isotropic (curves t) and 
anisotropic (curves a) scattering mechanisms are presented. 



III. RESULTS AND DISCUSSION 

The results of simulation, obtained for different 7 val- 
ues, are presented in Fig. ||. Let us consider, at first, the 
results for the isotropic scattering (dashed curves). They 
are in a good agreement with the results of numerical 
calculation carried out beyond the diffusion limitQ. This 
lends support to the validity of the method used. Here 
the results of calculation in the diffusion limit (6^1, 
7 <g; 1), obtained through the well-known expression 
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are shown too. As is clearly seen, the conditions, under 
which the expression (jj) works well, are very rigorous: 
even for 7 — 0.01 and 6 < 1 the formula (0) gives evi- 
dently larger value of Acr(6) than that in our simulation. 




FIG. 2. The dependences Aa{b) = Sa{b) - Sa{0) for the 
isotropic (dashed curves) and anisotropic (solid curves) scat- 
tering mechanisms, which are calculated with the angle de- 
pendences of scattering probabilities presented in the inset in 
Fig. h[_ Dotted curves are the results of diffusion approxima- 
tion 
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FIG. 3. The results of fitting of simulated curves (solid 
and dashed curves in Fig. hi) by expression (W). Dashed lines 
in upper figures correspond to -yfu = 7. 



The solid curves in Fig. || are the results of our sim- 
ulation for the anisotropic scattering mechanism. The 
closeness of these results to those, obtained for isotropic 
scattering, stands out. At first glane it seems to be 
surpising. The introduction of anisotropy in scattering 
process decreases the probability of returning to the start 
point area, because the scattering for small angles dom- 
inates in this case and in the average the particle moves 



to greater distance from the start point after each colh- 
sion. As is seen from (p , the decreasing of W has to lead 
to decreasing in 5a. However, this decreasing is compen- 
sated by the increasing of transport length I, which is 
averaged with the weight (1 — cos{9)) when integration 
runs over the scattering angle 9. Thus, the introduction 
of anisotropy of scattering does not dramatically change 
the negative magnetoresistance (at least up to 6 = 5). 

One can now attack the results of simulation as exper- 
imental ones. Let us describe our curves by expression 
(0) in the usual fashion. We have introduced a numerical 
multiplier (so-called prefactor) in (|4[) and used it and 7 
as fitting parameters. The results of such data process- 
ing are presented in Fig. ^. The curves labels show the 
magnetic field ranges, in which the fitting procedure has 
been carried out. In spite of strong difference between 
the simulation results and those given by (Q) (see dotted 
and other curves in Fig. 0), the fitting value of jfit is 
very close to 7, used in our simulation for both isotropic 
and anisotropic scattering mechanisms. The difference 
between 7 and jfa is less than 10 % for the isotropic 
scattering and 40 — 50 % for the anisotropic one. The 
value of prefactor is less than unity and decreases with 
increasing 7. 

Thus, (i) the use of expression (^ to determine the 
phase-breaking length (or time) from the magnetic field 
dependence of anomalous magnetoconductance in semi- 
conductor structures, where scattering is anisotropic, can 



give the error of about 40 — 50 %, (ii) the fact that the 
prefactor is less than 1 can be related not only with e- 
e interaction influence as it is frequently supposed, but 
with poor fulfilment of the condition I <C l,p and strong 
anisotropy of scattering as well. 
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